Search for: All records

Abundant epidemiological evidence links circadian rhythms to human health, from heart disease to neurodegeneration. Accurate determination of an individual’s circadian phase is critical for precision diagnostics and personalized timing of therapeutic interventions. To date, however, we still lack an assay for physiological time that is accurate, minimally burdensome to the patient, and readily generalizable to new data. Here, we present TimeMachine, an algorithm to predict the human circadian phase using gene expression in peripheral blood mononuclear cells from a single blood draw. Once trained on data from a single study, we validated the trained predictor against four independent datasets with distinct experimental protocols and assay platforms, demonstrating that it can be applied generalizably. Importantly, TimeMachine predicted circadian time with a median absolute error ranging from 1.65 to 2.7 h, regardless of systematic differences in experimental protocol and assay platform, without renormalizing the data or retraining the predictor. This feature enables it to be flexibly applied to both new samples and existing data without limitations on the transcriptomic profiling technology (microarray, RNAseq). We benchmark TimeMachine against competing approaches and identify the algorithmic features that contribute to its performance. 

Abstract An autonomous, environmentally-synchronizable circadian rhythm is a ubiquitous feature of life on Earth. In multicellular organisms, this rhythm is generated by a transcription–translation feedback loop present in nearly every cell that drives daily expression of thousands of genes in a tissue–dependent manner. Identifying the genes that are under circadian control can elucidate the mechanisms by which physiological processes are coordinated in multicellular organisms. Today, transcriptomic profiling at the single-cell level provides an unprecedented opportunity to understand the function of cell-level clocks. However, while many cycling detection algorithms have been developed to identify genes under circadian control in bulk transcriptomic data, it is not known how best to adapt these algorithms to single-cell RNAseq data. Here, we benchmark commonly used circadian detection methods on their reliability and efficiency when applied to single cell RNAseq data. Our results provide guidance on adapting existing cycling detection methods to the single-cell domain, and elucidate opportunities for more robust and efficient rhythm detection in single-cell data. We also propose a subsampling procedure combined with harmonic regression as an efficient, reliable strategy to detect circadian genes in the single–cell setting. 

The circadian clock is an evolutionarily-conserved molecular oscillator that enables species to anticipate rhythmic changes in their environment. At a molecular level, the core clock genes induce circadian oscillations in thousands of genes in a tissue–specific manner, orchestrating myriad biological processes. While previous studies have investigated how the core clock circuit responds to environmental perturbations such as temperature, the downstream effects of such perturbations on circadian regulation remain poorly understood. By analyzing bulk-RNA sequencing ofDrosophilafat bodies harvested from flies subjected to different environmental conditions, we demonstrate a highly condition-specific circadian transcriptome: genes are cycling in a temperature-specific manner, and the distributions of their phases also differ between the two conditions. Further employing a reference-based gene regulatory network (Reactome), we find evidence of increased gene-gene coordination at low temperatures and synchronization of rhythmic genes that are network neighbors. We report that the phase differences between cycling genes increase as a function of geodesic distance in the low temperature condition, suggesting increased coordination of cycling on the gene regulatory network. Our results suggest a potential mechanism whereby the circadian clock mediates the fly’s response to seasonal changes in temperature. 

Abstract Background During embryogenesis, the developmental potential of initially pluripotent cells becomes progressively restricted as they transit to lineage restricted states. The pluripotent cells of  Xenopus  blastula-stage embryos are an ideal system in which to study cell state transitions during developmental decision-making, as gene expression dynamics can be followed at high temporal resolution. Results Here we use transcriptomics to interrogate the process by which pluripotent cells transit to four different lineage-restricted states: neural progenitors, epidermis, endoderm and ventral mesoderm, providing quantitative insights into the dynamics of Waddington’s landscape. Our findings provide novel insights into why the neural progenitor state is the default lineage state for pluripotent cells and uncover novel components of lineage-specific gene regulation. These data reveal an unexpected overlap in the transcriptional responses to BMP4/7 and Activin signaling and provide mechanistic insight into how the timing of signaling inputs such as BMP are temporally controlled to ensure correct lineage decisions. Conclusions Together these analyses provide quantitative insights into the logic and dynamics of developmental decision making in early embryos. They also provide valuable lineage-specific time series data following the acquisition of specific lineage states during development. 

During embryogenesis, the developmental potential of initially pluripotent cells becomes progressively restricted as they transit to lineage restricted states. The pluripotent cells of Xenopus blastula-stage embryos are an ideal system in which to study cell state transitions during developmental decision-making, as gene expression dynamics can be followed at high temporal resolution. Here we use transcriptomics to interrogate the process by which pluripotent cells transit to four different lineage-restricted states: neural progenitors, epidermis, endoderm and ventral mesoderm, providing quantitative insights into the dynamics of Waddington’s landscape. Our findings shed light on why the neural progenitor state is the default lineage state for pluripotent cells, and uncover novel components of lineage-specific gene regulation. These data reveal an unexpected overlap in the transcriptional responses to BMP4/7 and activin signaling, and provide mechanistic insight into how the timing of signaling inputs such as BMP are temporally controlled to ensure correct lineage decisions. Together these analyses provide quantitative insights into the logic and dynamics of developmental decision making in early embryos. 

The most enigmatic of the canonical properties of circadian clocks is temperature compensation where circadian period length is stable across a wide temperature range despite the temperature dependence of most biochemical reactions. While the core mechanisms of circadian clocks have been well described, the molecular mechanisms of temperature compensation are poorly understood especially in animals. A major gap is the lack of temperature compensation mutants that do not themselves unambiguously affect the temperature dependence of the encoded protein. Here we show that null alleles of two genes encoding components of a complex important for translation of the core clock component period in circadian pacemaker neurons robustly alter the temperature dependence of circadian behavioral period length. These changes are accompanied by parallel temperature dependent changes in oscillations of the PER protein and are consistent with the model that these translation factors mediate the temperature-dependence of PER translation. Consistent with findings from modeling studies, we find that translation of the key negative feedback factor PER plays an instrumental role in temperature compensation. 

The univariate Kolmogorov-Smirnov (KS) test is a non-parametric statistical test designed to assess whether two samples come from the same underlying distribution. The versatility of the KS test has made it a cornerstone of statistical analysis across the scientific disciplines. However, the test proposed by Kolmogorov and Smirnov does not naturally extend to multidimensional distributions. Here, we present the fasano.franceschini.test package, an R implementation of the 2-D KS two-sample test as defined by Fasano and Franceschini (Fasano and Franceschini 1987) and provide multiple use cases across the scientific disciplines. The fasano.franceschini.test package provides three improvements over the current 2-D KS test on the Comprehensive R Archive Network (CRAN): (i) the Fasano and Franceschini test has been shown to run in O(n2) versus the Peacock implementation which runs in O(n3); (ii) the package implements a procedure for handling ties in the data; and (iii) the package implements a parallelized permutation procedure for improved significance testing. Ultimately, the fasano.franceschini.test package presents a robust statistical test for analyzing random samples defined in 2-dimensions. 

Ray et al . (Reports, 14 February 2020, p. 800) report apparent transcriptional circadian rhythms in mouse tissues lacking the core clock component BMAL1. To better understand these surprising results, we reanalyzed the associated data. We were unable to reproduce the original findings, nor could we identify reliably cycling genes. We conclude that there is insufficient evidence to support circadian transcriptional rhythms in the absence of Bmal1 . 

Circadian clocks play a key role in regulating a vast array of biological processes, with significant implications for human health. Accurate assessment of physiological time using transcriptional biomarkers found in human blood can significantly improve diagnosis of circadian disorders and optimize the delivery time of therapeutic treatments. To be useful, such a test must be accurate, minimally burdensome to the patient, and readily generalizable to new data. A major obstacle in development of gene expression biomarker tests is the diversity of measurement platforms and the inherent variability of the data, often resulting in predictors that perform well in the original datasets but cannot be universally applied to new samples collected in other settings. Here, we introduce TimeSignature, an algorithm that robustly infers circadian time from gene expression. We demonstrate its application in data from three independent studies using distinct microarrays and further validate it against a new set of samples profiled by RNA-sequencing. Our results show that TimeSignature is more accurate and efficient than competing methods, estimating circadian time to within 2 h for the majority of samples. Importantly, we demonstrate that once trained on data from a single study, the resulting predictor can be universally applied to yield highly accurate results in new data from other studies independent of differences in study population, patient protocol, or assay platform without renormalizing the data or retraining. This feature is unique among expression-based predictors and addresses a major challenge in the development of generalizable, clinically useful tests.